Rashba spin-splitting control at the surface of the topological insulator Bi2Se3 
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The electronic structure of Bi2Se3 is studied by angle-resolved photoemission and density func- 
tional theory. We show that the instability of the surface electronic properties, observed even in 
ultra-high-vacuum conditions, can be overcome via in-situ potassium deposition. In addition to ac- 
curately setting the carrier concentration, new Rashba-like spin-polarized states are induced, with 
a tunable, reversible, and highly stable spin splitting. Ab-initio slab calculations reveal that these 
Rashba state are derived from the 5QL quantum- well states. While the K-induced potential gradient 
enhances the spin splitting, this might be already present for pristine surfaces due to the symmetry 
breaking of the vacuum-solid interface. 

PACS numbers: 71.20.-b, 71.10.Pm, 73. 20. At, 73.22. Gk 



Topological insulators, with a gapless topological sur- 
face state (TSS) located in a large bulk bandgap, define 
a new quantum phase of matter [1-4 . Their uniqueness, 
and their strong application potential in quantum elec- 
tronic devices, stem from the TSS combination of spin 
polarization and protection from backscattering [5l |6]. 
Bi2Se3 is a three dimensional topological insulator, as 
theoretically proposed [7 and experimentally verified by 
angle-resolved photoemission spectroscopy (ARPES) and 
other surface sensitive techniques jSHTO]. Unfortunately, 
despite great effort in controlling the Se stoichiometry 
and with it the bulk carrier concentration [TJ, uninten- 
tional and uncontrolled doping seems to lead to a bulk 
conductivity that masks the surface electronic properties 
[12] . ARPES studies also have shown that cleaved sample 
surfaces and subsurfaces become progressively more elec- 
tron doped - even in ultra-high vacuum conditions - by 
either gas adsorption, or formation/migration of defects 
and vacancies [I3l|l4]. Lastly, the TSS might be com- 
pletely destroyed when exposed to air, hindering most 
attempts of material processing and characterization, as 
well as device fabrication. 

Developing new approaches to stabilize and control the 
surface of these systems is arguably the most critical step 
towards the exploitation of their topological properties. 
Some success has been obtained in inducing electron and 
hole surface doping by a combination of in-situ process- 
ing, such as material evaporation and radiation expo- 
sure [T5", T6 . The same TSS has also been fabricated 
on nanoribbons, which have large surface-to-volume ra- 
tio [17]. From a different perspective, carefully doped 
topological insulators can provide a platform to study 
the interplay between TSS and bulk electron dynamics, 
which has important implications for TSS control and 
exploring topological superconductivity [18j . 

In this paper, we present a systematic ARPES study of 



the evolution of the surface electronic structure of Bi2Se3 
as a function of time and in-situ potassium evaporation. 
The deposition of submonolayers of potassium allows us 
to stabilize the otherwise continually evolving surface 
carrier concentration. It also leads to a more uniform sur- 
face electronic structure, in which well-defined Rashba- 
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FIG. 1. (color online). (a,b) Time evolution of the ARPES 
dispersion of Bi2Se3 at 5xl0~^^ torr and T = 6K: (a) 3 hours 
after cleaving; (b) 34 hours after cleaving, (c) Exponential 
fit of the Dirac point (DP) binding energy versus time for 
6 and 300 K cleaves (both measured at 6K); saturation at 
AE^^ ^ 133 and 116 meV is reached in 46 and 22 hours, 
respectively. 
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FIG. 2. (color online). Evolution of the Bi2Se3 T—K electronic dispersion upon subsequent 0.5 min K-evaporation steps: (al-a8) 
ARPES image plots, (bl-b8) corresponding energy distribution curves (EDCs). The sample was kept at 5xl0~^^ torr and 6K. 



like states emerge from the continuum of parabolic-like 
states that characterizes the as-cleaved, disordered sur- 
faces. This approach provides a precise handle on the 
surface doping, and also allows tuning the spin splitting 
of the Rashba-like states. Our density functional theory 
(DFT) slab calculations reveal that the new spin-split 
states originate from the bulk-like quantum- well (QW) 
states of a 5-quintuple-layer (5QL) slab, as a consequence 
of the K-enhanced inversion symmetry breaking already 
present for the pristine surface of Bi2Se3. 

ARPES measurements were performed at UBC with 
21.2 eV linearly polarized photons on an ARPES spec- 
trometer equipped with a SPECS Phoibos 150 hemi- 
spherical analyzer and UVS300 monochromatized gas 
discharge lamp. Energy and angular resolution were 
set to lOmeV and ±0.1°. Bi2Se3 single crystals, grown 
from the melt (with carrier density n ~ 1.24 x 10^^ cm~^ 
pT] ) and by floating zone, were aligned by Laue diffrac- 
tion and cleaved and measured at pressures better than 
5xl0~^^torr and 6K, unless otherwise specified. No 
difference was observed for samples grown with differ- 
ent methods. Potassium was evaporated at 6K, with a 
6.2 A evaporation current for 30 second intervals [191 [20]. 
DFT calculations were performed using the linearized 
augmented-plane-wave method in the WIEN2K package 
[2T] , with structural parameters from Ref . 7 , We consid- 
ered stoichiometric slabs terminated by a Se layer on both 
sides, representing natural cleavage planes within this 
material. Spin-orbit interaction is included as a second 
variational step using scalar-relativistic eigenfunctions as 
a basis [21 ; exchange and correlation effects are treated 



within the generalized gradient approximation [22\ . 

The time evolution of the as-cleaved Bi2Se3 surface 
is shown in Fig[l] As typically observed by ARPES, 
and contrary to what is predicted by DFT for fully 
stoichiometric Bi2Se3 (Fig.|4|, even immediately after a 
6K cleave the Fermi level is not in the bulk gap; in- 
stead it crosses both TSS and parabolic continuum of 
bulk-like states. The pronounced time dependence of 
the data is exemplified by the variation of the Dirac 
point (DP) binding energy (AE^^), which increases 
from ~ 300 to 400 meV over 34 hours at 5xlO"^Horr 
and 6K [Fig.[l|c)]. An exponential fit of AE^^ versus 
time indicates that the saturation value E^^ 433 meV 
is reached 46 hours after cleaving. At variance with the 
time dependence of the TSS, the bottom of the parabolic 
continuum shifts down by only 30 meV in 34 hours, which 
provides evidence against the pure surface nature of the 
continuum. One should note that the pristine position of 
DP depends also on the cleave temperature: on a sam- 
ple cleaved at 300 K we found a 70meV deeper starting 
position for the DP, although the saturation value is ap- 
proximately the same as that of the 6K cleave [Fig.jljc)]. 

In our ARPES study, the surface time evolution re- 
sulted only in the deepening of Dirac cone (DC) and bulk 
continuum, as a consequence of the sample gaining elec- 
trons. Other effects, such as the reported appearance of 
a 2-dimensional electron gas (2DEG), were not observed 
^23j . More substantial changes are induced by the in-situ 
evaporation of potassium on the cleaved surfaces, also 
performed at 6 K to guarantee the highest stability. As a 
function of K-deposition time, three stages can be identi- 
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fied: Stage I- for moderate K deposition [up to 1 minute, 
Fig.|2|abl-ab3)], the DP moves to higher binding energy 
by electron doping and a sharper parabohc state appears 
at the edge of the bulk continuum, reminiscent of the 
proposed 2DEG [23 . Stage II: for intermediate K depo- 
sition [from 1 to 3 minutes, Fig.[2|ab4-ab7)], as the elec- 
tron doping further increases, the single parabola splits 
into a first pair - and then a second one also appears 
- of sharp parabolic states with an equal and opposite 
momentum shift away from the f point, as in a Rashba 
type [24 splitting [these states are labeled RBI and RB2 
in Fig. [3]^ a)]. Interestingly, the appearance of the sharp 
RBI and RB2 features is accompanied by a suppression 
of the bulk-like continuum. This emergence of a coherent 
quasiparticle dispersion from a continuum of incoherent 
spectral weight indicates that the evaporation of potas- 
sium leads to a progressively more uniform surface and 
subsurface structure. Stage III: for heavy K deposition 
[beyond 3 minutes, Fig.[2|ab8)], the bottom of RBI and 
RB2 as well as E^^ are not changing, indicating that the 
sample cannot be doped any further. The only notice- 
able effect is a small decrease of spin splitting for RBI 
and conversely an increase for RB2, perhaps stemming 
from a change in hybridization between the two Rashba 
pairs. As a last remark, during the entire K-deposition 
process the band velocity of the TSS close to the DP is 
3.2 ± 0.3 eVA, consistent with previous reports [25] . 

Before analyzing quantitatively the evolution of the 
various states upon K deposition, we address the question 
of the stability of this new surface versus time and tem- 
perature cycling. In Fig.[3|a-c) we compare the ARPES 
data from a 3 minute K-e vapor ated surface, as measured 
right after deposition and 30 hours later (during which 
the sample was kept at 6K). Other than a smaller than 
lOmeV shift of the bottom of RBI [Fig.[3];c)], ah spec- 
tral features including the TSS have remained exactly 
the same over the 30 hour interval. This is a remarkable 
stability, especially when compared to the 365 meV shift 
induced by the initial K deposition [Fig.|3|e)], and to the 
more than 100 meV shift observed versus time without 
any active surface processing [Fig.jljc)]. This approach 
might provide a new path to overcome the general insta- 
bility and self-doping problem of the surface of Bi2Se3, 
which represents one of the major shortcomings towards 
the fabrication of topological devices. Temperature ef- 
fects were studied by slowly warming up the sample, in 
which case K atoms diffuse and eventually leave the sur- 
face, reverting the material back to an earlier stage with 
lower K coverage. Indeed, as one can see by comparing 
Fig.jsjd) to Fig.[2|a4), a sample initially K evaporated 
for 4 minutes at 6K, and then measured at 220 K after 
a gradual 36 hour warming up, exhibits ARPES features 
similar to those obtained directly after a 1.5 minute K de- 
position at 6K. This implies that K deposition on Bi2Se3 
is also reversible, making it possible to fine tune surface 
doping, position of the DP, and Rashba spin splitting. 
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FIG. 3. (color online). ARPES T-K band dispersion from 
Bi2Se3 taken (a) immediately after a 3 min K evaporation, and 
(b) 30 hours later (the sample was kept at 5xl0~^^ torr and 
6 K the whole time) . As also emphasized by the comparison of 
the corresponding F point EDCs in (c), the evaporated surface 
is highly stable, (d) Band dispersion measured at 220 K after 
a slow 36 hour warming up on a sample initially K evaporated 
for 4 min at 6K; the comparison with the data in Fig.[2ja8) 
reveals the suppression of the K-induced carrier doping, (e- 
g) Evolution vs. K-evaporation time of: (e) binding energy 
variation for DP (AE^^) and bottom of RBI (A^^^^), as 
defined in (a); (f) sheet carrier density for DC (nf^) and RBI 
{^2D^)'i (g) variation of the DC Fermi wavevector {Akp^) 
and of the Rashba band splitting at Ep (Akp^-^). Empty 
symbols in (e-g) are for T = 6K and filled ones for T = 220K. 



We summarize in Fig.[3]^e-g) the K-evaporation evolu- 
tion of various parameters characterizing the T — K dis- 
persion of DC and Rashba states (empty symbols iden- 
tify 6K data, and the filled ones 220 K data). As evi- 
dent in Fig.js^e) from the variation of E^^ and bottom 
of RBI, the highest possible doping level is achieved ~3 
minutes into the K deposition, corresponding to AE^^ 
365 meV and A^^^^ :^ 150 meV. The K-induced change 
in surface electron density for the various states can be 
estimated from the relation n2D = ^fs /^bz ^uc be- 
tween the area of Fermi surface, Brillouin zone, and 
unit cell, without accounting for spin degeneracy given 
that all relevant states are spin split. Because at these 
electron fillings all FS's are hexagonal, this reduces to 
= ^F/2>/37r^, where hp is the Fermi wavevector along 
the F-^ direction of the BZ (as in Fig.|2|. After 3 
minute K evaporation the total sheet carrier density is 
77,2^^^3.64 X 10-^^cm~^ (0.162 electron/BZ), correspond- 
ing to the sum of the contributions from DC, and inner- 
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FIG. 4. (color online). DFT results for (a) /c^— projected bulk 
and (b) 5QL slab of Bi2Se3 (Ef is at energy 0). (c) Energy 
difference of the quantum-well (QW) slab states with respect 
to QWl, calculated for various thicknesses. The QW number 
increases with QL number, defining a bandwidth Wql that 
asymptotically approaches be fully /c^— projected bulk Wb = 
520 meV. The 3 minute K-deposition RB2-RB1 splitting of 
123±6meV is accurately reproduced by the 5QL results. 

and-outer RBI and RB2 (1.43, 1.60, and 0.61 x 10^^ cm-^ 
respectively). This value is to be compared to n*2D — 
3.87x10^2 cm-2 before K deposition (0.017 electron/BZ), 
which however only accounts for the DC, given the impos- 
sibility of estimating the contribution from the parabolic 
continuum. 

As a last point, from the Akp data presented in 
Fig.jsjg), and the dispersion of spin-split Rashba bands: 

EHkii) = Ep + ^±anh, (1) 

we can estimate the Rashba parameter aR = h^Ak^/2m'^ 
for RBI. The latter, which depends both on the value of 
spin-orbit coupling (SOC) and the gradient of the poten- 
tial dV/dz [26], reflects the size of the spin splitting in 
momentum space and is here controlled directly by the 
amount of K deposited on the as-cleaved surfaces. The 
largest RBI splitting is observed after 2.5-3 minute K 
evaporation and is anisotropic: Akp — 0.066A~^ along 
T-K, and O.OSOA"^ along f-M. Fitting the RBI disper- 
sion along T-K to Eq.^ we obtain Ef^ = 280±10meV, 
m* = 0.28±0.02me, and = 0.79±0.03eVA. This value 
is more than twice the Rashba splitting of the Au(lll) 
surface state (a^^ 0.33 eVA), and also larger than the 
one of the Bi(lll) surface state (a^^ :^0.56eVA) [27 . 

DFT calculations for bulk Bi2Se3, as well as slabs with 
varying number of QL's (Fig.|4|, provide a detailed ex- 
planation for our observations and some interesting in- 
sights. Each QL consists of 2 Bi and 3 Se layers alter- 
nating along the c axis, with one Se layer in the middle 
of the QL and the other two on either side. This forms 
a non-polar structure with a natural cleavage plane be- 
tween two adjacent Se layers belonging to different QL's. 
As shown in Fig.[4]^b) for the particular case of a 5QL 
slab, in addition to the TSS-DC there are 4 QW states. 



for a total of 5 states matching the number of QL's. As 
evidenced by the comparison with the fully /c^— projected 
bulk results in Fig.jlj^a), where the TSS is missing due to 
the absence of the surface, the slab QW states exhibit 
the same character and energy as the Bi-Se conduction 
band. However, they are discrete in nature due to quan- 
tum confinement, and span a narrower energy range than 
the corresponding bulk bandwidth Wb = 520 meV. The 
effective slab bandwidth Wql^ defined as the energy dif- 
ference between top and bottom QW states, is asymp- 
totically approaching the bulk Wb value [Fig.|4jc)]; for a 
proper correspondence with the bulk electronic structure 
a rather large number of QL's is needed (i.e., more than 
10 QL). Interestingly, the splitting between the DP and 
the different QW states is extremely sensitive to the num- 
ber of QL's. For 5QL we obtain 346 meV QWl-DP and 
126 meV QW2-QW1 splittings [Fig.[4](c)], which closely 
match the 3 minute K-deposition values 380 ± 50 meV 
and 123±6meV for RBl-DP and R_B2-RB1, respectively 
[as defined from the EDC's at the f point in Fig.jsja)]. 

This analysis leads to several important conclusions: 
(i) The RBI and RB2 states that emerge from the 
parabolic continuum are of the same Bi and Se char- 
acter than those obtained, in the same energy range, in 
bulk Bi2Se3 calculations. However, because of the ob- 
served lack of kz dispersion [23 , and the almost perfect 
match with QWl and QW2 energy positions for 5QL 
[as also emphasized by the comparison between calcula- 
tion and 3 minute K-deposition datapoint in Fig.[4]^c)], 
these states should be more appropriately thought as 
the quantum-confined analog of those bulk states, asso- 
ciated with a band-bending over a 5QL subsurface region 
(47.7 A). While this region is somewhat disordered - ei- 
ther in its depth and/or carrier concentration - on the 
as-cleaved surfaces, it becomes progressively better de- 
fined upon K evaporation, (ii) Potassium, in addition to 
doping carriers, also induces a change in dV/dz^ which in 
turn provides a very direct control knob on both band- 
bending depth and spin splitting of the Rashba states. 
(in) In light of the extent of the subsurface band-bending 
region, these quantum-confined states should affect more 
than just surface sensitive experiments. For instance, 
Rashba spin-split states might have to be accounted for 
in the interpretation of transport data even from pristine 
surfaces, although with a much smaller splitting induced 
solely by the symmetry breaking vacuum-solid interface. 
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